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The Earth is subjected to a uniform flux of very-high-energy (VHE, E > lOOGeV) cosmic rays 
unless they are obscured by an object, such as the Moon, in which case a deficit or shadow is 
created. Since cosmic rays are charged this deficit is deflected by the Earth’s magnetic held, 
enabling the rigidity of the obstructed cosmic rays to be determined. Measurement of the relative 
deficits of different species enables the positron fraction and the antiproton ratio to be measured. 
The April 15, 2014 lunar eclipse was visible with the VERITAS Cherenkov telescopes, which 
allowed (with special UV bandpass Alters) 74 minutes of direct observations of the Moon and 
the associated deficit in the cosmic-ray flux. The results of this observation are presented. In 
addition VERITAS has been conducting a series of observations by pointing close to a partially 
illuminated Moon, with a reduced photomultiplier tube high voltage and UV bandpass Alters. We 
present the technique developed for these observations and their current status. 
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1. Inroduction 

Since their discovery over 100 years ago, cosmic rays have been the subject of intense re¬ 
search, yet their origin remains a mystery. Key evidence to help solve this problem will come from 
the energy spectrum of cosmic rays and the relative proportions of the different cosmic ray types, in 
particular the ratio of antimatter to matter. In this work we are looking at two important measures: 
the positron fraction (c+/(c+ + e^)) in the range ~ 0.5 to ~ ITeV and the antiproton ratio (p/p) 
in the ~ 1.5 to ~ 3TeV range. This is done using the Earth-Moon ion spectrometer (EMIS), pio¬ 
neered for use with lACTs by the ARTEMIS collaboration [1] and also being used by the MAGIC 
collaboration to measure the positron fraction [2, 3]. Observations are subject to constraints on 
Moon elevation and illumination, thus observing time is limited to a few hours per year, this time 
includes observations during a Eunar eclipse on 2014-04-15. 

2. Positron Fraction & Antiproton Ratio 

Since the first conclusive result of an increase in the positron fraction from the PAMEEA 
collaboration [4] and especially with the improved statistics and energy range provided by AMS- 
02 [5] ,the positron fraction has been of particular interest. This rise was largely unexpected and 
deviates markedly from “traditional” models as shown in Eigure 1 where a typical GALPROP [6] 
result is overlaid. The key question has been what will happen at higher energies where a measure 
of the fraction can greatly constrain the existing models. 



Figure 1: Current measurements of the positron fraction from AMS-02 [5], PAMEEA [4] and the 
Fermi-CPW [7]. Overlaid is a typical GAEPROP prediction [6] of the positron fraction showing the 
divergence above ~10GeV. Data is from [8]. 

The antiproton ratio is well measured up to a few hundred GeV (with data from AMS-02 
recently reducing the statistical uncertainty and increasing the energy range). However, above 
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1 TeV there are only upper limits, as can be seen in Figure 2. At present, none of these limits are 
constraining on typical models that fit the existing antiproton ratio measurements at lower energies. 
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Figure 2: The antiproton ratio (from [9]) showing the upper limits from a variety of experiments 
above ITeV, the solid lines are predictions based upon pure secondary production [10] and the 
dot-dash line includes contributions from additional sources of p at cosmic ray accelerators [11]. 


3. The Earth-Moon Ion Spectrometer 

The Earth-Moon ion spectrometer (EMIS) is a technique to determine the relative fluxes of 
cosmic rays of different rigidities {R = p/q where p is the particle momentum, q is the charge). The 
basis of the technique is that the cosmic ray flux is uniform across the sky. The only exception is if 
there is an object obstructing this flux, creating a deficit or shadow, such as the Moon. Since cosmic 
rays are deflected by the Earth’s magnetic field, with the amount of deflection depending upon the 
rigidity of the cosmic ray, this deficit will also be deflected. When viewing this shadow from the 
ground (Eigure 3) it creates a band that passes across the Moon, with one charge deflected one way, 
one the other and with lower rigidity particles deflected further than higher rigidity particles. The 
ratio between the strength of the positive and negative shadows can be used to determine the ratio 
between positive and negative species. By applying analysis cuts, hadronic and leptonic showers 
can be differentiated, allowing for separate ratios to be determined for the two. 

4. Feasibility of Observing with VERITAS 
4.1 Bright Moonlight Observations 

In order to conduct these observations the telescopes must be pointed close to a partially il¬ 
luminated Moon. The Moon is a significant light source that will increase the background light, 
reducing the signal-to-noise ratio, increasing the trigger rate (and thus the deadtime) and poten¬ 
tially damaging the PMTs. To reduce its impact two techniques were employed, the high voltage 
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Figure 3: A diagram showing the cosmic ray shadow the amount of deflection (shown for various 
rigidities as dotted circles) depends on the position in the sky. For a falling Moon in the west, the 
positive shadow is at a lower elevation. 

supplied to the PMTs was reduced (RHV) and Schott glass UG-11 filters (UVF) were employed 
that only pass UV light, reducing the amount of light entering the PMTs. These methods and their 
impact on the performance of the array are discussed in detail elsewhere in this conference [12]. 
For the normal Moon shadow observations, a combined approach with both RHV and UVF is used. 
During the 2014 lunar eclipse, due to the reduced lunar brightness we were able to observe using 
only UVF. 

4.2 Shadow Position 

In order to conduct the observations we need to know where the shadow is. The deflection 
of a particle is a function of the magnetic field, with the amount of deflection depending upon the 
rigidity of the particle. In order to find this position we used the code TJ2010 [13]. The shadow axis 
is defined using the Moon and the position of the ±750GV shadow. Observations are conducted 
with the centre of the camera pointed ±0.5° perpendicular to this line. Since the shadow is the 
width of the Moon (0.5°) and with an angular resolution of 0.1° for leptonic showers, this pointing 
means the shadow is fully contained in one half of the camera leaving the other half for background 
estimation. 

4.3 Crab Nebula Observations 

On 2013-03-19 the Crab Nebula was ~3° from a 45% illuminated Moon, this provided an 
opportunity to test the sensitivity of the array in the proposed observing conditions (UVF -i- RHV). 
Observations were conducted for all of the time the Crab Nebula was above 40° elevation, with a 
total live time of 107.5 minutes. These data show that even with both RHV and UVF the VERITAS 
array is still sensitive to /-rays, though at a much reduced level of 6.82 a/V hours using the stan¬ 
dard analysis cuts. For comparison, normal observations have a sensitivity of around 30 a/ y/hours 
at these elevations. 
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4.4 Observation Simulation 

In order to test the observation’s feasibility and to develop the analysis techniques, a Monte- 
Carlo model was generated assuming an energy threshold of 500 GeV for leptons and 1.5TeV for 
hadrons. This simulation (smoothed with a Gaussian kernel of 0.1°) is shown in Figure 4, the total 
number of counts falling within on field-of-view (FoV) is ~2.2 x 10^. This corresponds to 30 hours 
of observations at a reconstruction rate of 10 Hz. The proton shadow is clearly visible cutting off 
sharply due to the power law spectrum just before ITV shadow position, along with the fainter 
helium shadow which extends through the ITV shadow position to the 500GV position due to 
its factor or two lower rigidity. The angular resolution (which is much larger at 0.4° for hadronic 
showers) extends the shadow in all directions and will be a major constraint in this analysis since it 
will reduce the sensitivity by smearing out the signal and reducing the area available for estimating 
the background. 
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Figure 4: Moon shadow simulation, the number of events within the FoV (large dashed circle) is 
approximately equivalent to 30 hours of observations. The shadow shape is influenced by the large 
angular resolution (0.4°), high energy threshold (1.5TeV) and the source spectrum (a power law of 
index -2.7) which results in the shadow lying close to the Moon in a “blob” rather than an band. 


If selection cuts are applied, the survival fraction is 0.3 for leptons and 0.001 for hadrons, a 
simulation with the counts reduced to these fractions is shown in Figure 5. The electron deficit 
just above the threshold energy is visible at a rigidity of -500 to —700GV, though a lot weaker 
than the proton deficit in the previous simulation. With the lower statistics, statistical fluctuations 
produce comparable deficits elsewhere in the simulation, though the one associated with the elec¬ 
tron shadow is present in all simulations, thus the 30 hours of observations can be considered a 
minimum. The positron deficit is not visible and requires significantly more time to detect using 
this simple analysis. However, neither simulation has made use of the reconstructed energy to aid 
in the analysis, the addition of this parameter will significantly improve the sensitivity. 
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Figure 5: A simulation of the expected skymap after gamma/hadron selection cuts are applied. 
The number of events within the FoV is approximately equivalent to 30 hours of observations. 

4.5 Feasibility Study Conclusions 

This study shows that it may be feasible to detect the shadow with VERITAS. Taking tens 
of hours to detect either shadow (potentially hundreds for the positron shadow) it will require a 
multiyear investment as a maximum of 15 — 20 hours of observing time can be expected in any one 
year (less with bad weather). However, since an upper limit on the positron fraction in the 0.5 to 
1 TeV range could potentially provide useful constraints, a detection of the positron shadow may 
not be required. Observations and support work are ongoing. 

5. Normal EMIS Observations 

Observations taken during the 2013-14 observing season totalled 326 minutes (out of a maxi¬ 
mum possible of ~900 minutes). These results do not show the presence of the cosmic ray shadows 
on their own (and would require significant more time to do so) but they have been used to verify 
that the tracking and the analysis method is working correctly. 

6. Lunar Eclipse Observations 

On 2014-4-15 a total lunar eclipse was visible from the VERITAS facility. 74 minutes of 
observations were taken pointing directly at the Moon in the elevation range 45 to 50°. During 
the majority of the time the currents were low (less than background stars that do not normally 
cause an effect on the analysis). The last 14 minutes had slightly higher currents, reaching levels 
comparable with bright stars towards the end of the data taking and thus this time has not been used 
in this analysis. Quality cuts were also chosen to minimise the impact of any optical brightness (at 
the cost of statistics and a raised energy threshold) by using a large Size (number of digital counts 
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in the cleaned image) cut and requiring that a minimum of three telescopes are used in the image 
reconstruction. Up to the event reconstruction stage the standard analysis packages are used, after 
that custom methods are employed (with an independent cross check) to plot the positions of the 
showers in a coordinate system centred on the Moon, but with the axes defined by the midline of the 
shadow (Parallel) and its perpendicular (Perpendicular). The Parallel distance is then normalised 
so that the distance between the 1 TV shadow position and the Moon is constant. The results are 
shown in Figure 6 where a deficit due to the shadowing of the Moon is clearly visible slightly offset 
from the Moon along the Parallel axis. 



Parallel [deg] 


Figure 6: The 20725 showers that passed the quality cuts are plotted relative to the Moon in the 
coordinate system described in the text with a 0.1° Gaussian smoothing applied. The deficit due to 
the Moon is slightly offset from the center of the camera. 

To test for deflection, a fit was conducted using a 2D Gaussian over the central part of the FoV. 
This fit gives a centroid location of (Parallel, Perpendicular) = (0.09°, -0.06°), a = (0.19°, 0.20°). 
It should be noted that the energy threshold is estimated at above 4TeV at which energy the shift 
is expected to be 0.2°. Further work is ongoing to determine whether there is any evidence for 
deflection. If selection cuts are applied to examine the electron shadow then only 32 counts remain 
and no conclusions can be drawn. 

7. Conclusion 

VERITAS is using the EMIS technique to measure the positron fraction in the energy range 
0.5 TeV to ITeV; a measurement which has been of great interest to the community over the last 
few years. This measurement will require many tens of hours, (if not hundreds for the positron 
shadow) taking many years to accumulate the data. Observations conducted during a lunar eclipse 
and pointed directly at the Moon show a clear deficit in the cosmic ray flux due fo fhe shadowing 
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of the Moon. Work is ongoing to determine whether there is any evidenee for a defleetion in the 
shadow position from the Moon. 
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